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Abstract 

O ' We show that the critical mass M c = 8tt of bacterial populations in two dimensions in 

the chemotactic problem is the counterpart of the critical temperature T c = GMm/Aks 
of self-gravitating Brownian particles in two-dimensional gravity. We obtain these critical 
O I values by using the Virial theorem or by considering stationary solutions of the Keller-Segel 

model and Smoluchowski-Poisson system. We also consider the case of one dimensional 
systems and develop the connection with the Burgers equation. Finally we discuss the 
evolution of the system as a function of M or T in bounded and unbounded domains 
in dimensions d = 1, 2 and 3 and show the specificities of each dimension. This paper 
\& ■ aims to point out the numerous analogies between bacterial populations, self-gravitating 

! Brownian particles and, occasionally, two-dimensional vortices. 

\6 
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1 Introduction 

In many fields of physics, astrophysics and biology, one is confronted with the description of the 
evolution of a system of particles which self-consistently attract each other over large distances. 
One difficulty and richness of the problem arises from the long-range nature of the potential 
of interaction [TJ. This is the case, for example, in biology in relation with the process of 
chemotaxis [2j . The chemotactic aggregation of bacterial populations (like Escherichia coli) or 
amoebae (like Dictyostelium discoideum) is usually studied in the framework of the Keller-Segel 
model [3] which describes the collective motion of organisms that are attracted by a chemical 
substance (pheromone) that they produce themselves. The Keller-Segel (KS) model involves a 
drift-diffusion equation describing the evolution of the concentration of the bacteria in the gra- 
dient of concentration of the secreted chemical. In the simplest formulation, the concentration 
of the chemical is related to the concentration of bacteria by a Poisson equation (this is valid in 
a limit of large diffusivity of the chemical and for sufficiently large concentrations) pE]. These 
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equations have been studied by applied mathematicians who obtained rigorous results for the 
existence and unicity of the solutions and for the conditions of blow-up, modeling chemotactic 
collapse, in different dimensions of space [5]. In particular, in d = 2, there exists a critical mass 
M c (independent on the size of the domain) above which the system collapses and forms a Dirac 
peak. Alternatively, for M < M c , the system spreads to infinity in an unbounded domain or 
tends to a stationary state in a bounded domain. 

Gravity is another example of long-range attractive potential of interaction. In a series of 
papers, Chavanis & Sire [61 El El El [TOj EH [121 E31 Ej have studied a model of self-gravitating 
Brownian particles in various dimensions of space. In statistical mechanics, this model is 
associated with the canonical ensemble in which the temperature is fixed. A lot of analytical 
results have been obtained and an almost complete description of the system, for all the phases 
of the dynamics (pre-collapse and post-collapse), has been given in the overdamped limit of 
the model. In that limit, the evolution of the density of the self-gravitating Brownian gas is 
governed by the Smoluchowski-Poisson (SP) system. The Smoluchowski equation is a drift- 
diffusion equation of a Fokker-Planck type. For self-gravitating particles, the gravitational 
potential inducing the drift is produced by the density of particles through the Newton- Poisson 
equation. It turns out that the Smoluchowski-Poisson system is isomorphic to the simplified 
version of the Keller-Segel model of chemotaxis provided that the parameters are suitably 
re-interpreted, as discussed in [11]. In particular, for the 2D self-gravitating Brownian gas, 
there exists a critical temperature T c (independent on the size of the domain) below which 
the system collapses and forms a Dirac peak. This is the counterpart of the critical mass of 
bacterial populations. For T > T c , the system evaporates in an unbounded domain or tends to 
a statistical equilibrium state in a bounded domain. 

The object of this paper is to emphasize the parallel between these two systems. In Sees. 
[2] and El we use the Virial theorem to derive the critical mass of bacterial populations and the 
critical temperature of self-gravitating Brownian particles in two dimensions. In Sec. HI we show 
that these critical values can also be obtained by considering stationary solutions of the Keller- 
Segel model and Smoluchowski-Poisson system. In Sec. we consider the one dimensional 
problem and point out the connection between the Smoluchowski-Poisson system (or the Keller- 
Segel model) and the Burgers equation. We use this analogy to provide the general solution of 
these equations in d — 1 in bounded and unbounded domains. Finally, in Sec. El we provide 
a summary of the results obtained by Chavanis & Sire [HI El El El QUI Ell E21 E31 Ell f° r se ^" 
gravitating Brownian particles and adapt them to the context of chemotaxis to clearly show 
the link between these two problems. 



2 The Keller-Segel model 

The dynamical evolution of biological populations like bacteria, amoebae, cells... that are 
attracted by a substance that they emit themselves, is often described by the Keller-Segel (KS) 
model |3J. In its simplest form, it can be written as 

(1) ^ = J DAp- X V.(pVc), 



(2) Ac = -Xp. 

Equation ([T]) is a drift-diffusion equation for the cell density p(r,t). The diffusion term takes 
into account the erratic motion of the cells (like in Brownian theory) and the drift term with 
X > takes into account the chemotactic attraction (one could also consider the case \ < 
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where the secreted substance is a noxious substance, like a poison, so that chemotaxis is 
repulsive) . It is directed along the gradient of concentration c(r, t) of the secreted chemical. 
In the simplest formulation [I], the production of the chemical by the cells is described by 
a Poisson equation (jSJ). This is valid in a limit of high diffusivity of the chemical and for 
sufficiently large concentrations (see Appendix The Keller-Segel model ([I])-© must be 
supplemented by appropriate boundary conditions. A first physical boundary condition is that 
the current J = D Vp — XP^c is parallel to the boundary of the domain (or vanishes at infinity 
in an unbounded domain) so that the total mass is conserved. On the other hand, when dealing 
with the reduced Keller-Segel model (JT])-([2]), one usually assumes that c(r, t) is the solution of 
the Poisson equation (j2J) in an infinite domain with the usual gauge condition. This yields 
c ( r >*) = (rf=fe / P( r V)l r - r'l'^W (in d ^ 2) or c(r,t) = ±J p(r',t)ln|r - r'\dr' (in 
d = 2) whenever the domain containing the bacteria is finite (box) or infinite (here Sd denotes 
the surface of a unit sphere in (i-dimensions) . A physical discussion of the boundary conditions 
and of the limitations of the reduced Keller-Segel model ([I])-© is provided in Appendix lAl 

We now proceed in deriving an exact relation that is similar to the Virial theorem in astro- 
physics (see Sec. [3]). As discussed in [13J, it is convenient to take the origin of the system of 
coordinates at the center of mass which is a fixed quantity. Multiplying Eq. ([T]) by XiXj and 
integrating over the entire domain, we get 

Introducing the "tensor of inertia" 

(4) = J pXiXj dr, 
and integrating the second term by parts twice, we find that 

(5) 1 -^L = DMS^ + X W t3 , 

where M = J pdr is the total mass of cells and Wij = Wji is the "potential energy tensor" [13J: 

(6) W t3 = J px^dr. 

If the system is confined within a box, we have to account for boundary terms. The first 
integration by parts in Eq. ([3]) yields a residual term 



dS, 



(7) / V • M DV P - XP Vc)] ir = f XiXl (DV P - XP Vc) ■ 

where dS is the surface element normal to the frontier of the confining box. By virtue of the 
conservation of mass, the diffusion current in Eq. (TjQ) is always perpendicular to the surface 
vector and consequently the term ([7j) vanishes. The second integration by parts yields 

(8) - J Dp(xj5 k i + Xi5 kj ) dr = - j> Dp(xj5 ki + Xi5 kj ) dS k . 

Therefore, the general form of the Virial theorem for the chemotactic problem, taking into 
account boundary terms, is 

(9) \lt = Dm ? + xWl > ~\j Dp ( X i dSi + Xi dS ^- 
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By contracting the indices, we get the scalar Virial theorem 

1 dl f 

(10) -— = dDM + xW ii -Dj>pr-dS, 

where 

(11) I = J pr 2 dr, 
is the moment of inertia and 

(12) W u = [ pr-Vcdv, 



is the Virial. If the density p& is uniform on the edge of the box (this is the case at least for a 
spherically symmetric system), we get 

(13) j) p r • dS = pb j) r • dS = pi J V • r dr = dp b V, 
where V is the volume of the confining box. Thus, 

(14) ~ = dDM + X W U - dD Pb V. 

On the other hand, adapting the results of the Appendix of Ref [13] to the present situation, 
we find for d ^ 2 that Wu — (d — 2)W where W = — | J pcdr is the "potential energy" . In that 
case, we get 

(15) ~^ = dDM+{d-2) X W -D J> pr ■ dS, (d^2). 
2 dt j 

For d = 2, we find instead that Wu = — AM 2 /(47r) and we obtain 

(16) I^ = 2DM-^^-D Ipr-dS (d = 2). 
At equilibrium (1 = 0), the Virial theorem fflUj) reduces to 

(17) dDM + xWu = D j pr-dS. 
For d = 2, we have 

(18) 2DM- = D j pr ■ dS. 

In usual situations, the term in the r.h.s. is positive (this is at least the case for an axisymmetric 
distribution of particles in a disk where it is equal to 2Dp b irR 2 ). Since p > 0, a sufficient 
condition is that r • dS > on each point of the boundary. This criterion is independent on 
the distribution itself and only depends on the domain shape. When <f> p r • dS > [J, the above 
relation implies that a necessary condition for the existence of steady solutions is that 

19 M < M c = — — . 



lr There are probably situations where this inequality is violated, in which case we expect that the value of 
the critical mass will differ from Eq. (|19p . In such situations, boundary effects should play a prominent role. 
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If we introduce dimensionless parameters (or take D — x — A = 1), the critical mass is simply 
M c = 8tt. In terms of the critical mass (|T9|) . we can rewrite the Virial theorem (|T6j) as 

(20) l^ = ^L( Mc -M)-Dj pr-dS (d = 2). 

For M > M c , we have / < e < so that the moment of inertia goes to zero in a finite time. 
This implies that the system collapses to a Dirac peak at r = (the center of mass) in a finite 
time. In an unbounded domain, assuming that p decreases more rapidly than r~ 2 for r — > +oo, 
the foregoing relation reduces to 

(21) 5§ = ^<*-"> 

This relation shows that, in an unbounded domain, there can be stationary solutions (J = 0) 
only when M = M c . In that case, the moment of inertia is conserved. For M < M c , the 
moment of inertia diverges with time I(t) — > +oo so that the system evaporates. For M > M c 
the moment of inertia goes to zero /(£)—► in a finite time so that the system forms a Dirac 
peak at r = in a finite time. Since the case M = M c lies at the frontiere between these two 
regimes (evaporation or collapse) we expect that, when M = M c , the system evolves either 
towards a Dirac peak (collapse) or a completely spread profile (evaporation). 

The existence of a critical mass for chemotaxis in d = 2 has been found by various authors 
[igUgimHamiSIlHSllSU, using different arguments. 



3 The Smoluchowski-Poisson system 



In the mean field approximation and in a strong friction limit, the dynamics of a one-component 
self-gravitating Brownian gas is described by the Smoluchowski-Poisson (SP) system 



(22) 



dt 



1 fk B T 



m 



- Vp + pV$ 



(23) A$ = S d Gp. 

The Smoluchowski-Poisson system has been studied by Chavanis & Sire [SI El 13 UUI [III 
IT2"| [TBI EH] i n different dimensions of space (see also extensions in [2"2"l |2"5] for more general 
equations of state p = p(p))- Since the temperature is fixed, the relevant statistical ensemble 
is the canonical ensemble. The boundary conditions are: (i) the current J = ^^Vp + pV$ is 
parallel to the boundary in a finite domain and vanishes at infinity in an unbounded domain 
(so that the total mass is conserved), (ii) The gravitational potential is given by $(r, t) = 
-r^r J p(r',t)\r - r'|-( d - 2 W (in d ^ 2) or $(r,t) = —G J p(r', t) In |r - r'\dr' (in d = 2), 
using the usual Gauge condition. 

It is clear at first sight that the Smoluchowski-Poisson system (I22p - (l2"3"j) is isomorphic to the 
simplified Keller-Segel model (JI])-(j25 provided that we make the correspondances 

(24) D = ^, x =- c=-$, X = S d G. 

Furthermore, the structure of the solutions depends on a single dimensionless parameter [7] 
which can be written: 

_ pGMm _ X X M 
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We note, in particular, that the concentration of the chemical c(r, t) in the chemotactic problem 
is the counterpart of the gravitational potential $(r, t) for self-gravitating Brownian particles 
(with the opposite sign). Due to this analogy, the results derived for the SP system can be 
applied to the KS model and vice versa. However, due to the different notations (and also 
because chemotaxis and gravity are studied by different communities) this connection is not 
always made. It can be therefore of interest to put the two models in parallel, as we do here. 

The Virial theorem for the SP system has been derived in [13] (see also Appendix [B]). It 
can be written 

(26) -A = dNk B T+(d-2)W-dPV, (d ^ 2) 

2 at 

where W = | J pQdr is the potential energy and we have defined 

(27) P = ±-j pr . d s^ 

where p(r) = p(r)kBT/m is the local pressure of an isothermal gas. In the case where p = pb is 
constant on the boundary, we have P = pb- In dimension d = 2, the Virial theorem takes the 
simple form [IB] : 

(28) ^§ = 2Nk B T - ^f- - 2PV, (d = 2). 
At equilibrium (1 = 0), we obtain 

CM 2 

(29) 2Nk B T — = 2PV, 

Since P > in usual circumstances (this is at least the case for an axisymmetric system in 
a disk where P = pbk B T/m), the above relation implies that a necessary condition for the 
existence of steady solutions is that 

nn nn GMm 

30 T>T C = ——. 

Ak B 

If we introduce dimensionless parameters (or take G = M = m = k B = £ = 1), the critical 
temperature is simply T c = 1/4. In terms of the critical temperature (130]) . the Virial theorem 
can be rewritten 



(31) ^ = 2Nk B (T-T c )-2PV. 

For T < T c , the moment of inertia goes to zero in a finite time so that the system collapses to 
a Dirac peak at r = in a finite time. At equilibrium (1 = 0), we obtain 

(32) PV = Nk B (T-T c ). 

This is the equation of state of the two-dimensional self-gravitating gas at statistical equilibrium 
in the thermodynamic limit N — ► +oo with rj = GMm/k B T fixed (in that limit, the mean field 
approximation is exact [21]). This equation of state (and its extension to finite N systems) has 
been obtained by various authors using different methods [25], [26] [27] US] [HI | 
In an unbounded domain, the Virial theorem ( |3TI) reduces to 

(33) ^f t =2Nk B (T-T c ). 



6 



This relation shows that, in an unbounded domain, there can be stationary solutions (/ = 0) 
only at the critical temperature T = T c . In that case, the moment of inertia is conserved. For 
T > T c , the moment of inertia diverges I(t) — > +00 so that the system evaporates. For T <T C 
the moment of inertia goes to zero I(t) — > in a finite time so that the system forms a Dirac 
peak at r = in a finite time. Defining the mean-squared radius of the cluster through the 
relation (r 2 ) = I/M, we can rewrite Eq. (1331) after integration in the form (r 2 ) = 4D (T)t + (r 2 ) 
where 

(34) jD (T) = M(l-T c /T), 

is an effective diffusion coefficient. For T>T C when gravitational effects become negligible, the 
Smoluchowski equation ( 122]) reduces to a pure diffusion equation and the diffusion coefficient is 
given by the Einstein formula D(+oo) = ksT j^m. However, Eq. ( 1341) shows that the diffusion 
is less and less effective as temperature decreases and gravitational effects come into play. In 
particular, the effective diffusion coefficient becomes negative for T <T C indicating finite time 
collapse. 

The existence of a critical temperature for self-gravitating systems in d = 2 has been found 
by various authors, in different contexts, using different arguments [2§1 EDI ED E21 El EE]- A 
similar critical temperature (whose value is negative) appears in the statistical mechanics of 
point vortices in two-dimensional hydrodynamics [331 EH ES] . 

4 The equilibrium density profile in d = 2 

4.1 Chemotaxis of bacterial populations 

The stationary solutions of the KS model ([I])-© are such that 

(35) p = Aeo c . 

This is similar to the Boltzmann distribution in statistical mechanics provided that we interprete 
T e ff = D/x as an effective temperature and — c(r) as a potential. This distribution can be 
obtained by extremizing the functional 

If D f 

(36) F[p] — — - J pcdr + —J plnpdr, 

at fixed mass. This functional is the Lyapunov functional of the KS model. It satisfies F < 
and F — if, and only if, the density is given by Eq. (1351) . It can also be interpreted as an 
effective free energy F = E — T e ffS associated with the Boltzmann entropy [36]. A steady state 
of the KS model is an extremum of F at fixed mass. Furthermore, it is linearly dynamically 
stable if and only if it is a (local) minimum of this functional [37J. The equilibrium state is 
obtained by substituting Eq. (135]) into the Poisson equation ([2} yielding 

(37) Ac = -XAeo c . 

This is similar to the Boltzmann-Poisson equation appearing in astrophysics (in the statistical 
mechanics of stellar systems [27, 38J and in isothermal models of stars [39]) and in vortex 
dynamics (in the statistical mechanics of 2D point vortices [35]). If we restrict ourselves to 
axisymmetric solutions, this equation can be solved analytically in d = 2. The 2D axisymmetric 
Boltzmann-Poisson equation (157]) is solved in Appendix O Here we use a slightly different 
method working directly on the mass profile. 
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Considering axisymmetric density profiles, introducing the accumulated mass M(r, t) = 
Jq p(r' , t)2irr' dr' , using the relation J*; = — AJ ^ 7 |^' t ' > (equivalent to the Gauss theorem) and in- 
troducing the variable u = r 2 , it is found that the KS model ([I])-© is equivalent to the single 
partial differential equation 

DM 4jDm <9 2 M x X M dM 
dt du 2 7r du 

The stationary profiles satisfy 

(39) uM" + ^-MM' = 0. 

Using mM" = («¥')' - M', 2MM' = (M 2 )' and M(0) = 0, we obtain after integration the first 
order differential equation 

(40) uM' = M ( 1 - -^M 



Since M' > 0, a necessary condition for the existence of steady states is that 
(41) M<M C 8nD 



Solving Eq. fj4TJ]) when this condition is fulfilled, we find that the equilibrium mass profile is 
given by 

Kr 2 

(42) M(r) 



1 i Xr 2 ' 
1 M c 



where is a constant of integration determined by the boundary conditions. 

(a) In a disk of radius R, the constant K is determined by the condition M(R) = M. This 
yields 

M (-) 2 

(43) M(r) - 



1 M I I M /r\2' 



The corresponding density profile is 

( 44 ) p(0 = r^T^ 



1 + M (^) 2 
1 ~ M C -M\R> 



2 ■ 



Typical density profiles are plotted in Fig. [T]for different values of the ratio M/M c . The relation 
between the mass and the central density is 

(AV s M c M/M c 

(45) Po = ^JT^- 

It is plotted in Fig. [2j The central density increases as M — > M c . For M = M c , the density 
profile is a Dirac peak: p(r) = M c S(r). 

(b) In an unbounded domain, we see from Eq. (l4"2"j) that the condition M{r) — * M for 
r — > +oo requires that M = M c . Therefore, there exists steady state solutions only for M = M c . 
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MM =0.99999 

c 




Figure 1: Density profile (normalized by M c /irR 2 ) as a function of r/R for different values of 
the total mass M. For M — > M c , the density profile tends to a Dirac peak. 




1 1 1 1 1 1 1 1 1 1 1 1 1 

-2 2 4 6 8 10 

ln(P ) 

Figure 2: Relation between the mass M/M c and the central density po (normalized by M c /ttR 2 ) 
in the chemotactic problem. There exists a unique solution for each value of M < M c . Since 
there is no turning point in the series of equilibria M(p ), these steady solutions are stable 
(minima of free energy F at fixed mass M) [7]. 
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They form an infinite family of solutions parameterized by the constant K or equivalently by 
the central density po — K/tt. Their density profile is 



(46) P(r) = - P ° „ N2 . 



1 + 



wpor 2 
M c 



It is suggested in Sec. |2] that, among all the distributions of this family, the only stable 
stationary solution of the KS model with M = M c is the Dirac peak corresponding to po — > +00. 
The dynamical formation of this Dirac peak is shown analytically in [13J. Another possible 
evolution is an evaporation. 



4.2 Self-gravitating Brownian particles 

We give, for comparison, the equivalent relations for the gravitational problem. The stationary 
solutions of the Smoluchowski equation fl22|) are described by the Boltzmann distribution 



(47) p = Ae-f 3 ™*. 

This is equivalent to the condition of hydrostatic equilibrium Vp + pV$ = for an isothermal 
equation of state p(r) = p{r)k,BT /m [13J. The Lyapunov functional of the SP system is the 
Boltzmann free energy 

(48) F[p] = - [ p$ dr + k B T [ — In — dr. 

2 J J m m 

The steady states are obtained by solving the Boltzmann-Poisson equation 

(49) A$ = S d GAe-^. 

They are linearly dynamically stable if, and only if, they are minima of free energy at fixed mass. 
This is consistent with a condition of thermodynamical stability in the canonical ensemble. In 
d = 2, the SP system is equivalent to a single partial differential equation for the mass profile 

dM d 2 M dM 

(a) In a bounded domain, there exists steady states if, and only, if 

(KU T^T GMm 

51 T>T C = — . 

Ak B 

The mass profile is given by 



(52) M{r) 



M (i) 2 



1 _ Ic X , _^_(2L) 
± T j. -r T _ Tc K R ) 

and the corresponding density profile by 

M 1 



( 53 ) P( r ) = — - , 7 r 

1 + ^M-V 

1 T T _ Tc \ R l 



T 



The relation between the temperature and the central density is 
(54) po 1 M 



vri? 2 l-f 
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At T = T c , the density is a Dirac peak p(r) = M5(r). 

(b) In an unbounded domain, there exists steady solutions only for T = T c . They are 
parameterized by the central density po — K/tt. The mass and density profiles are 

(55) M(r) = Kr l - , p(r) = ? . 

! + f^ 2 (1 + ^ 2 ) 

These solutions have the same value of the free energy F = —NkBT c [l +ln(7r/iV)] independent 
on po- On the other hand, their moment of inertia is infinite (except the solution with po = +00 
for which 1 = 0). Since / is conserved at T = T c , they cannot be reached from a generic initial 
condition with < I < +00. It is suggested in Sec. [3] that the only stable stationary solution 
of the SP system with T = T c is the Dirac peak corresponding to p — > +00. The dynamical 
formation of this Dirac peak is shown analytically in [T3]. The system ejects a tiny amount 
of mass at large distances in order to accomodate for the conservation of /. Another possible 
evolution at T = T c is an evaporation. 

In the above formulae, we can pass from biological notations to gravitational notations by 
using the relation 



(56) 



T M, 



c 



T r . M 



resulting from the correspondances (12^1) . Therefore, the mass M in chemotaxis plays the role of 
the inverse temperature (3 = l/fegT in gravity. We note also that, in 2D gravity, the condition 
( 15"Tj) for the existence of a steady solution above a critical temperature T c (for a given mass 
M), can be converted into a condition on the mass (for a given temperature T) 

(57) M < M c = 

making the link with the biological problem even closer. 



5 Analogy with the Burgers equation in d = 1 
5.1 The case of an infinite domain 

In d = 1, it is shown in [7] that the Smoluchowski-Poisson system is equivalent to a single 
differential equation 

(») ^ = ^ + GM™ 

at m ox z ox 

for the integrated density M(x, t) = 2 J* p(x', t)dx' . The density and the gravitational potential 
are related to the mass profile by d x M = 2p(x,t) and d x Q = GM(x,t). The stationary states 
of Eq. (158j) can be obtained analytically (see Sec. 15.21) . In an infinite domain, they are given 
by 

(59) M(x) = Mtanh(x/H), p(x) 



cosh (x/H) 

where the length scale and the central density are determined as a function of the temperature 
by 

. , 2k B T M GM 2 m 

{60) = 'GMm 1 P0 = 2H = ^f- 
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There exists a unique steady state solution for each value of the temperature. Introducing 
rescaled parameters, or taking equivalently £ = ks = m = G = M = 1, the dynamical equation 
(1581) takes the form 



(61) 

Setting v(x,t) = —M(x,t) and v 
equation 

(62) 



dM 



8 2 M ,BM 
T-^ + M- 



dt dx 2 dx 
- T, it becomes equivalent to the one dimensional Burgers 



dv dv 
dt dx 



d 2 v 
dx 2 



which has been studied in connection with hydrodynamical turbulence jlQ] and cosmology [H] . 
In this analogy, the mass plays the role of the velocity (v = —M), the density the role of the 
velocity increment (v' = — 2p) and the temperature the role of the viscosity [y = T). The 
steady solution of Eq. (1621) is 



(63) 



v(x) = — tanh(x/2i/). 



In ID turbulence, this solution describes a single shock [ID] while in the context of self- 
gravitating Brownian particles it describes the equilibrium mass profile fl59l) of a one dimensional 
isothermal gas [13J. 

The collapse dynamics of the Smoluchowski-Poisson system at T = has been studied in 
[7. 8j where an analytical solution describing the formation of a Dirac peak has been obtained. 
The collapse is self-similar and a Dirac peak is formed in the post-collapse regime [7J [8]. This 
Dirac peak is precisely the limiting form of the stationary solution (l59l) - (!60l) at T = 0. We note 
that the solution obtained in [TIE] also describes the formation of a singular shock in inviscid ID 
turbulence (y = 0) since the equation for the mass profile of self-gravitating Brownian particles 
in d — 1 is isomorphic to the Burgers equation. For v ^ 0, the Burgers equation ( |62l) has a 
well-known explicit solution. As observed by Hopf [12] and Cole [13] , the change of variables 
v = —d x ip and ip = 2z/ln# transforms the nonlinear Burgers equation into the heat equation 
thereby leading to the explicit solution 



(64) 



ip(x, t) = 2v In 



+oo 



exp 



1 

27 



^o(g) - 



(x — q)' 
2T~ 



dq 



where ipo(x) = ip(x,0). Returning to the notations of the initial problem, we note that if) 
represents the gravitational potential $. This provides the general solution of the SP system 
in d = 1 dimension in an unbounded domain for T^O. Acedo [2] has constructed an explicit 
analytical solution of the ID Smoluchowski-Poisson system from the general formula If 
we take the limit T — > and use steepest descent technics, we obtain from Eq. (|64j) : 



(65) 



if>(x,t) = sup 



^o(g) - 



[X 



qf 



2t 



(T-0). 



A general method to solve this equation is described in [JT] in the cosmological context. Inter- 
estingly, these results can also have applications for the ID Smoluchowski-Poisson system. 



5.2 The case of a finite domain 

According to Eq. (J58"j) . the stationary solution of the SP system in d — 1 satisfies the differential 
equation 



(66) 



M" + ^MM' = 0. 
k B T 
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x/R 

Figure 3: Density profile (normalized by M/2R) as a function of x/R for different values of the 
temperature rj = /3GMmR. For T — > 0, the density profile tends to a Dirac peak. 



Integrating this equation once and using M(0) = and M'(0) = 2po to determine the constant 
of integration, we find that 



(67) 

This can be rewritten 
(68) 



, ,, Gm , 

M + — M = 2p 

2k B T H0 



a 2 - M 2 2fc s T 
This is easily integrated in 



Gm 



(69) 



M(rc) = atanh 



/ Gma 

\2k B r 



p{x) 



Po 



cosh 



2 / Gma 
2k B T 



X 



In an infinite domain, using the fact that M(x) — > M for x — > +oo, we find that a = M and 
we obtain the relations (|59l) -( l60l) . In a bounded domain, using the fact that M(R) = M we 
find that a is solution of 



(70) 



M = a tanh 



/ GmaR 



This equation implicitly determines the central density po m terms of the temperature T. These 
results can also be obtained by solving the ID Boltzmann-Poisson equation (see Appendix iDl) 
which is equivalent to Eq. (166]) [7]. The equilibrium density profile for different values of the 
temperature is represented in Fig. [3] and the relation between the temperature and the central 
density is represented in Fig. HI We have used the notations defined in Appendix [Dl 

Introducing rescaled variables or taking equivalently £ = k B = m . = G = M = R ', = 1, 
the dynamical equation for the mass profile in d = 1 is given by Eq. (16 ip with the boundary 
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Figure 4: Relation between the normalized inverse temperature rj = (3GMmR and the central 
density po (normalized by M/2R). There exists a unique solution for each value of the temper- 
ature. Since there is no turning point in the series of equilibria i](po), these steady solutions 
are stable (minima of free energy F at fixed mass M) [TJ. 

conditions M(0,t) = 0, M(l,t) = 1 and M(x, 0) = M (x). With the change of variables 
v(x,t) = —M(x,t), v = —d x ip and ip = 2Tln# of Sec. 15.11 we find that the function 9(x,t) 
satisfies the diffusion equation 



(71) 

with the boundary conditions 



d6 d 2 9 
— = T . 

dt dx 2 



(72) 



9'(0,t) = 0, 9'(l,t) = ^9(l,t), 9 (x) = e &J7"o(v)* 



2T 



Using the method of separation of variables, we find that the general solution of these equations 



is 



(73) 



9(x, t) = a e TA «* cosh(A x) + ^ a n e TA "' cos(A r 



n=l 



where Ao is the solution of tanh(Ao) = l/(2TAo) and X n are the solutions of tan(A r 
— l/(2TA n ). The coefficients are determined from the initial condition 



(74) 



9q(x) = a cosh(Aoa;) + a n cos(A n a;). 



n=l 



Using the relations satisfied by Ao and A n , the following identities can be obtained by standard 
calculations 



(75) 



1 / -i % M 1 + 4T 2 A2-2T X 

cos(A„x) cos{X m x)dx = o n 

2(1 + 4T 2 A 2 ) 
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"1 

(76) / cosh(A x) cos(X n x)dx = 0, 

o 



(77) /'cosh^Aox)^- 4 ^ 1 + 2T 



2(4T 2 A 2 -1) • 

Then, we find that the coefficients a n are given in terms of the initial condition by 

(78) a = ^^3__}L J 6o{x) cos h(\ x)dx, 

2(1 + 4T 2 A 2 ) f 1 

(79) a n = /rTnn - = — / 9n(x) cos(\ n x)dx. 

The foregoing equations provide the general solution of the SP system in <i = 1 in a bounded 
domain. For t — > +oo, we find from Eq. (1731 that M(x,t) — > M(x) = 2TX tanh(Aox) with 
1 = 2TAotanh(Ao) which returns the stationary solution ( J69l) -(|70l). 

5.3 Generalization to other dimensions 

Note, finally, that the Smoluchowski-Poisson system in d dimensions is equivalent to the single 
differential equation [7]: 

, . dM k B T fd 2 M d-ldM\ G dM 



dt m \ dr 2 r dr J r d 1 <9r ' 

for the mass profile M(r, t) = J Q r p(r' ,t)Sdr' d ~ l dr' '. With the change of variables x = r d /d, the 
foregoing equation becomes 

(81) * dM - keT (dx) 2( ^ 1)d2M I CM 8M 

dt m dx 2 dx 

Introducing rescaled variables, or using equivalently £ = kB = m = G = M = R=l, and 
setting u (x, t) = —M(x,t), this can be put in the form of a Burgers equation: 

/„„s 9v dv , . <9 2 v 

(82) ¥ + "fc = J/(x) a?' 

with a position dependent viscosity 

(83) v[x) = T(d x) 2 -^. 

The viscosity is a pure power- law scaling like 1, x and x 4//3 in d — 1,2 and 3 dimensions 
respectively. For T = 0, we obtain the inviscid Burgers equation in any dimension of space. 

In conclusion, it is interesting to note that the spherically symmetric Smoluchowski-Poisson 
system is connected to the ID Burgers equation which arises in many domains of physics. This 
increases the interest of the self-gravitating Brownian gas model [B]. We note however that, 
in our case, v(x,t) = —M(x,t) is a negative quantity for x > while in ID turbulence the 
velocity v can take positive and negative values. 
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6 Dynamics of bacterial populations and self-gravitating 
Brownian particles 

The Virial theorem (I26p can be used to obtain general results on the dynamics of self-gravitating 
Brownian particles without solving the equations of motion. However, the explicit resolution of 
the Smoluchowski-Poisson system, as considered by Chavanis & Sire [6],[7l[8|[^[T0 t [TT |[T2l[T3| [TI]. 
gives more precise information on the evolution of the density profile of the particles. In this 
section, we provide a summary of these results in the case of self-gravitating Brownian particles 
and transpose them to the chemotactic problem by using the notations of biology. This should 
reinforce the link between these two topics that are studied by different communities. 

6.1 Self-gravitating Brownian particles 

First, consider a two-dimensional self-gravitating Brownian gas enclosed within a circular box 
of radius R. (a) Regular steady states exist for T > T c and their density profile is known 
analytically, see Eq. (J53l) . They are global minima of the free energy ()48l) at fixed mass. Using 
the free energy as a Lyapunov functional, we deduce that the SP system relaxes towards these 
steady states. The particles are confined by the box since P > (i.e. p(R) > 0). (b) For 
T = T c , the steady state is a Dirac peak p(r) = M5(r) containing the whole mass so that 
P = (i.e. p{R) = 0). The dynamics of the SP system at T = T c has been studied in Sec. 
IV.B. of |7j. The collapse is self-similar with a scaling function similar to Eq. (153]) . and the 
Dirac peak containing the whole mass is formed for t — > +oo. The central density increases 
exponentially rapidly as p (t) ~ e^ 1 . (c) For T < T c , there is no steady state and the system 
undergoes gravitational collapse. The dynamics of the SP system at T < T c has been studied 
in Sec. IV. C. of [TJ. For T < T c the collapse is not exactly self-similar. In a finite time t co u, 
the system develops a Dirac peak containing a fraction T/T c of the total mass M, surrounded 
by a halo whose tail decreases as p ~ r~ a ^ with a(t) converging extremely slowly to a = 2 for 
t — > t co u, like 2 — a(t) = a/2 In In p (t ) / In Po (t) with po(t) ~ (t co u — £) _1 , so that the evolution 
looks self-similar. For t — > t co u, the density profile behaves like p(r,t) (T/T c )MS(r) + p(r,t) 
where the evolution of p(r,t) is studied in [7j. From the Virial theorem Eq. (13TT) . we note that 
/ < e < for T < T c so that the moment of inertia tends to its minimum value / = in a finite 
time t en d. This corresponds to the formation of a Dirac peak at r = containing the whole 
mass M. Since this final state is different from the structure obtained at t = t co u, this means 
that the evolution continues in the post-collapse regime t co u < t < t end . In this post-collapse 
regime, the Dirac peak formed at t = t co u accretes the mass of the surrounding halo until all 
the mass is at r = at t en d- This regime has not yet been fully described. 

Consider now a two-dimensional self-gravitating Brownian gas in an infinite domain, (a) For 
T > T c , the particles have a diffusive motion (evaporation) with an effective diffusion coefficient 
given by Eq. (134H increasing linearly with the distance T — T c to the critical temperature. This 
evaporation process is self-similar and it has been analytically studied in Sees. IV.B. and IV. C. 
of [13] for T>T C and T — ► T 6 + respectively, (b) For T = T c , the effective diffusion coefficient 
vanishes D(T C ) = 0, so that the moment of inertia is conserved: 1 = 0. The dynamics of the 
SP system at T = T c has been studied in Sec. D. of [13] . There exists an infinite family of 
steady state solutions (151)1) parameterized by the central density p but they have / = +oo (for 
Po < +oo) and cannot be reached by the system for generic initial conditions (with / < +oo). 
An analytical dynamical solution has been found in which the system forms a Dirac peak of 
mass M — e for t — ► +oo and ejects a tiny amount of mass e < 1 at large distances so as 
to satisfy the moment of inertia constraint (if I ^ initially). Note that in an unbounded 
domain [13], the central density diverges logarithmically with time, po{t) ~ lnt, while the 
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divergence is exponential in a bounded domain [7J. Another possible evolution at T = T c is an 
evaporation, (c) For T < T c , the effective diffusion coefficient in Eq. (jMl) is negative, implying 
finite time blow-up. In particular, (r 2 ) = for t en d = (r 2 ) /4:\D(T)\ where (r 2 ) is calculated 
at t = from the center of mass. This leads to the formation of a Dirac peak p(r) = M5(r) at 
r = containing the whole mass. The solution obtained in [Tjj in a bounded domain probably 
describes the collapse of the core accurately but, in the absence of a confining box, the collapse 
is accompanied by an unlimited expansion of the halo. In the post-collapse regime, the Virial 
theorem indicates that all the matter falls at r = in a finite time. 

For comparison, let us describe the situation in d — 3. When the system is confined within 
a spherical box of radius R, there exists a critical temperature fcsT c = f 52R which depends on 
the box radius R (contrary to the case d = 2). (a) For T > T c , there exists equilibrium states. 
They have a density contrast p(0)/p(R) < 32.1. These equilibrium states are metastable^, i.e. 
local minima of the free energy (1481) at fixed mass [BJ[7J. For T = T c , the equilibrium state is not 
a Dirac peak (contrary to the case d = 2). The density profile is regular with a density contrast 
p(0)/p(R) = 32.1. (b) For < T < T c , there is no equilibrium state and the system undergoes 
a self-similar collapse p(r,t) = po(t) f (r / r (t)) leading to a finite time singularity at t — t co u. 
The collapse time behaves with the distance to the critical point as t co u =t*{j] — ri c )~ l l 2 where 
1] = j3GMm/R and £* = 0.91767702... [9]. This pre-collapse can be described analytically in 
d > 2 and the invariant density profile is known exactly; the density decreases as r~ 2 at large 
distances P,[7j. The central density increases with time as p ~ {t co u —t)~ l and the core radius 
decreases as r (t) ~ (t co u — t) 1 / 2 . Therefore, this pre-collapse regime does not create a core 
since the mass M (t) ~ p r d at r = goes to zero as M (t) ~ {t co u — t)^' 2 ^ 2 at t = t co u. 
However, a Dirac peak is formed in the post-collapse regime for t > t co u [8]. The mass of the 
Dirac peak increases as M (t) ~ (t — ico«) f° r t — > t^ oll and the residual density obeys a 

backward dynamical scaling p(r,t) = po(t)g(r/r (t)) with p ~ {t — tcoii) 1 and r ~ {t — t co u) l l 2 . 
The Dirac peak accrets all the mass in an infinite time t —>■ +oo. (c) For T = 0, the system 
undergoes a self-similar collapse and develops a finite time singularity at t — t co u which can be 
studied analytically in any dimension. The invariant density profile is solution of an implicit 
equation and the density decreases as r~ 2d ^ d+2 ' > at large distances [HIE]- The central density 
increases with time as p ~ {poUi~ an d the core radius decreases as r Q (i) ~ (t co u — t)( d + 2 )/ M . 
Therefore, this pre-collapse regime does not create a core since the mass M (t) ~ {tcoii ~ t) d ^ 2 
at r = goes to zero at t co u. However, a Dirac peak is formed in the post-collapse regime 
for t > tcou [H [8] . The mass of the Dirac peak increases as 

and the residual density obeys a backward dynamical scaling p(r,t) = Po(t)g(r/ro(t)) with 
Po ~ {t — tcou)^ 1 and r ~ (t — tcott) 2d ■ At T = 0, the Dirac peak accretes all the mass in 
a finite time t en d = 1/d [5]. 

For 3D self-gravitating Brownian particles in an unbounded domain, there is no steady 
state. The system can either collapse (as in a bounded domain) or evaporate. The evaporation 
process has been treated in Sec. V. of [13J and the collapse in [61 UJ. The choice between 
collapse or evaporation depends on a complicated notion of basin of attraction that is function 
of the initial condition. 

Finally, in d = 1 the problem can be mapped on the Burgers equation which has the explicit 
solution (16~41) in an unbounded domain and (1731) in a bounded domain. 

2 Since there is no global minimum of free energy for an isothermal self-gravitating gas, these equilibrium 
states may not always be reached by the system. For some peculiar initial conditions, the system may rather 
collapse and form a singularity (Dirac peak) with infinite free energy as in the case T < T c (see below). This 
depends on a complicated notion of basin of attraction as illustrated in [6] . However, for generic initial conditions 
with T > T c , the SP system usually converges towards a "gaseous" steady state with smooth density profile 
(local minimum of free energy). 
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6.2 Chemotactic aggregation of bacterial populations 

The reduced Keller-Segel model (EE])-© and the Smoluchowski-Poisson system (I22"I) - (I23"1) are 
isomorphic so we can directly transpose the results obtained for self-gravitating Brownian par- 
ticles to the biological context. By reformulating these results with notations appropriate to 
chemotaxis, we aim to facilitate the comparison with the numerous results obtained in chemo- 
taxis by applied mathematicians [5j. Since we are just transposing the results of the previous 
section to a different context, we shall give less details in their description. 

First, consider the two-dimensional axisymmetric Keller-Segel model ([I])-© in a circular 
domain of radius R. (a) Regular steady states exist for M < M c and their density profile 
is known analytically, see Eq. (|44p . (b) For M = M c , the steady state is a Dirac peak 
p(r) = M c 5(r) containing the whole mass. The collapse is self-similar and the Dirac peak 
is formed for t — > +00. The central density increases exponentially rapidly as po(t) ~ e^*. 
(c) For M > M c , there is no steady state and the system undergoes chemotactic collapse. For 
M > M c the collapse is not exactly self-similar. In a finite time t co u, the system develops a Dirac 
peak of mass M c surrounded by a halo whose tail decreases as p ~ r - "^ with a(t) converging 
extremely slowly to a = 2 for t — > t co u, so that the evolution looks self-similar. This corresponds 
to p(r, t) — > M c S(r) + p(r, t). This non self-similar collapse has been studied in detail in [TJ. 
Some results had been obtained earlier by Herrero & Velazquez [18] but they are different from 
those of [7] . The reason is not well-understood but it may be that the solutions constructed in 
[18] . while mathematically correct, are unstable. By contrast, the solution of [7] shows a very 
good agreement with direct numerical simulations of the reduced Keller-Segel model in d = 2 
(as explained previously, the work [7] is presented for self-gravitating Brownian particles but 
the results also apply to biological populations). From the Virial theorem Eq. (1201) . we note 
that / < e < for M > M c so that the moment of inertia tends to its minimum value 1 = 
in a finite time t en d- This corresponds to the formation of a Dirac peak at r = containing 
the whole mass M. Therefore, we expect a post-collapse regime for t co u < t < t en d leading 
ultimately to a Dirac peak p(r) = M5(r) containing all the mass at t en d- To our knowledge, 
this post-collapse regime has not yet been fully characterized. 

Consider now the case of 2D systems in an infinite domain, (a) For M < M c , I(t) — > +00 
for t — ► +00 and the system evaporates. This regime has been studied in [13] . (b) For M = M c , 
dl/dt = so that the moment of inertia is conserved. It is shown analytically in [13] that the 
system forms a Dirac peak for t — ► +00 and ejects a tiny amount of mass at large distances 
so as to satisfy the moment of inertia constraint. In an unbounded domain [13], the central 
density diverges logarithmically with time, po(t) ~ hit, while the divergence is exponentially 
fast in a bounded domain [TJ. For M = M c , the system can also evaporate, (c) For M > M c , 
I(t) = in a finite time t en d so the system collapses to a Dirac peak containing all the mass 
in a finite time. The evolution is expected to be similar to the case of a box-confined system 
studied in [7J. 

In d = 3, in a bounded domain, there exists a critical mass M c = 31.7^ which depends on 
the box radius R (contrary to the case d = 2). This is the equivalent of the Emden temperature 
ksT c = 252R m astrophysics, (a) For M < M c , there exists metastable equilibrium states, 
i.e. local minima of the functional f l36l) . to which the system is attracted for generic initial 
conditions, (b) For M > M c , there is no equilibrium state and the system undergoes a self- 
similar collapse leading to a finite time singularity at t = t co u. The collapse time behaves 
with the distance to the critical point as t co u = ^(77 — r] c )^ 1 ^ 2 where 77 = XMx/^DR and 
U = 0.91767702... [9J. This pre-collapse regime can be described analytically [HI [7J. The 
density profile decreases as r~ 2 . The central density increases with time as po ~ (tcou — t)^ 1 
and the core radius decreases as r (t) ~ (t co u—t) 1 / 2 . A Dirac peak is formed in the post-collapse 



18 



regime for t > t co u. The mass of the Dirac peak increases as M (t) ~ (t — t co u) l l 2 . The Dirac 
peak accrets all the mass for t — > +00 |8J. (c) The peculiar case D = (no diffusion) where 
the particles are driven solely by the chemotactic drift is equivalent to the case T = for 
self-gravitating Brownian particles described in Sec. 16.11 

In an unbounded domain in d = 3, there is no steady state. The system can either collapse 
(as in a bounded domain) or evaporate. The evaporation process has been treated in [T5] 
and the collapse in [6], [7] . The choice between collapse or evaporation probably depends on a 
complicated notion of basin of attraction that is function of the initial condition. 

Finally, in d — 1 the problem can be mapped on the Burgers equation which has the explicit 
solution f l64l) in an unbounded domain and ( 1731) in a bounded domain. 



In this paper, we have stressed the analogy between the chemotaxis of bacterial populations 
and the dynamics of self-gravitating Brownian particles. In particular, in d — 2 dimensions, 
we have shown that the critical mass of bacterial populations M c is the counterpart of the 
critical temperature T c of self-gravitating brownian particles. These analogies are not well- 
known because these topics (chemotaxis and gravity) are usually studied by very different 
communities and the self-gravitating Brownian gas model has been introduced only recently 
in physics [6]. Yet, we think that the inter-relation between these disciplines is important to 
develop and the present paper is a step in that direction. In particular, we have obtained the 
value of the critical mass by using a relation which turns out to be equivalent to the Virial 
theorem in astrophysics. We think that many other connections can be made between the two 
disciplines, and this will be considered in future works. 

We have also qualitatively discussed the dynamical evolution of a self-gravitating Brownian 
gas (or a chemotactic system) in different dimensions of space by presenting a synthesis of the 
results obtained in Chavanis & Sire [6l [7J [8j [9j [TOl [TTJ [121 EES E] • This gives a clear picture of 
the collapse dynamics of a spherically symmetric system to a single cluster. If we come back to 
the general problem which does not need to be spherically symmetric, we expect that several 
collapses will take place at different locations of the domain (if sufficiently large). Each collapse 
will be described by the spherical solution that we have found. But the resulting clusters will 
themselves have a non-trivial dynamics and will "merge" together so that their number will 
decrease with time until a single Dirac peak containing the whole mass remains at the end. 
In that case, the evolution toward the final Dirac peak is progressive. This problem shares 
some analogies with the dynamics of vortices in 2D decaying turbulence jlSJ H6] although the 
equations of motion are of course different. The analogy with 2D turbulence may be interesting 
to develop. The aggregation of clumps in our gravitational Brownian model [6J could also 
be studied by exploiting the analogy with the Burgers equation [JT]. These are directions of 
investigation that we plan to explore in the future. 

A A remark on the boundary conditions 

When we consider the more general Keller-Segel model [3j HI [5j [11] 



7 Conclusion 



(84) 





— = D'Ac + ap- be, 
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the boundary conditions are the Neumann conditions 

(86) Vp • n = Vc • n = 0, 

where n is a unit vector normal to the boundary of the box. Equation (I85p describes the 
evolution of the concentration of the chemical. The chemical diffuses with a diffusion coefficient 
D', is created by the bacteria at a rate a and is degraded at a rate —b. With the Neumann 
boundary conditions, there is no current of particles (bacteria and chemical) accross the box 
so we have p = c = outside the box. From Eqs. and (IBS]) , we note that the average 

concentration of bacteria p(t) = M/V = p is conserved. From Eqs. fl85|) and (IBS]) , we note 
that the average concentration of chemical satisfies dc/dt + be = ap so that c(t) = (ap /6)(l — 
e~ bt ) + c e~ bt tending to c(+oo) = ap /b at equilibrium. Following Jager & Luckhaus [Ij, we 
set a = XD' and consider the limit D' — > +oo with A ~ 1. This leads to 

(87) ^ = jDA p- x V(pVc), 



(88) Ac=-X(p-p ). 

This model is still well-posed mathematically with the Neumann boundary conditions (IBS]) . 
We emphasize that only the gradient of concentration Vc enters in Eqs. ( 186]) - (IBB]) . Therefore, 
the concentration c itself is un-determined when we make the above-mentioned approximations 
since it is obtained only within an additive constant. When the density blows up (chemotactic 
collapse) so that p(r,t) 3> p , it is justified to consider the model ([I])-© where Eq. ( IBB] ) 
is replaced by a Poisson equation. It is only in that case (large diffusivity of the chemical 
D' — > +oo and high concentration p(r, t) ^> p of the bacteria) that the Keller-Segel model 
for the chemotaxis becomes equivalent to the Smoluchowski-Poisson system for self-gravitating 
Brownian particles. Some authors have studied the model ([I])-© on general grounds, i.e. not 
necessarily being an approximation valid when p(r, t) ^> p . We stress, however, that this 
model (0Q)-(I2]) is not well-posed mathematically with the Neumann boundary conditions (IBS"]) . 
Indeed, integrating Eq. (T5]) and using the divergence theorem, we have § Vcn dS = —AM ^ 0, 
so that the Neumann boundary conditions (IBS"]) cannot be satisfied in that case. One way to 
circumvent this difficulty is to use the boundary conditions defined in Sec. |2j We must keep 
in mind, however, that these boundary conditions do not determine the physical concentration 
Cphys{ r ,t) but only a "field" c(r,t) that has the same gradient (this field can take positive or 
negative values). Furthermore, this field does not satisfy the requirement Vc • n = so it is 
not clear how close it is related to the solution of the problem (IBS]) - (IBB]) . Finally, although the 
problem (JTJ)- (J2J) is well-posed mathematically with the boundary conditions of Sec. [2] we now 
find that c(r, t) ^ outside the box. This is a physical problem because there is no reason 
why the chemical should exit the material box (unless it has porous properties). Of course, the 
practical solution is to consider the concentration c(r, t) inside the box only, and ignore that 
field outside. Another possibility is to impose c = on the boundary of the domain (Dirichlet) 
and c(r, t) = outside. In the case of self-gravitating Brownian particles, the concentration 
— c(r,t) is replaced by the gravitational potential $(r,t). We can have $(r, t) ^ outside 
the box enclosing the particles because the gravitational field "traverses" the box (a material 
body enclosed within a container creates a gravitational force outside this container) so that 
the boundary conditions defined in Sec. [3] are natural. We conclude therefore that the reduced 
chemotactic model CQ)-(j2]) is not very well-posed physically compared with the initial problem 
( IB4]) - (IBS]) . By contrast, the Smoluchowski-Poisson system (T22l - (f23j) is rigorous on a physical 
point of view as the gravitational potential is always solution of the Newton-Poisson equation 
( 123]) with the boundary conditions given in Sec. [31 
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B Overdamped Virial theorem for spherically symmet- 
ric systems 

In this Appendix, we show that the overdamped Virial theorem (1261) can be established very 
simply in the case of spherically symmetric systems. In that case, the moment of inertia can 
be expressed as a function of the mass profile in the form 



(89) 



f R dM f R 
I(t) = -jr(r, t) r 2 dr = MR 2 - 2 / M(r, t)r dr, 
Jo v r Jo 



where the second equality follows from an integration by parts. Taking the time derivative of 
Eq. ( l89i) . inserting Eq. ( IHUl) and using simple integrations by parts, we obtain 

This can be rewritten in the form 

where we have used P = ksT 'p(R) jm and V = ^S^-R^. Noting finally (see the Appendix of 
that 



f R GM(r,t)dM t . , 

( 92 ) Wu = ~l -^=^^ t)dr ^ 

we obtain Eq. (1261) . In particular, for d = 2 we immediately deduce Eq. (1281) from Eq. (19TT) . 

C Explicit solution of the 2D Boltzmann-Poisson equa- 
tion 

In this Appendix, we explicitly solve the Boltzmann-Poisson equation (|37|) in d = 2 for axisym- 
metric solutions. The density (|35|) can be rewritten 

(93) p = p e^ c - c °), 

where p an d Cq are the values of the density at the center of the domain. The Boltzmann- 
Poisson equation takes the form 

(94) l± (rpi = -\p ei^\ 

r dr \ dr J 

Introducing ip = — -^(c — c ) and £ = (XxPo/ Dy/ 2 r, we obtain 



(96) ijj(0) = ^'(0) = 0. 
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With the change of variables t = In £ and ip — 2 In £ — z, Eq. (|95|) can be rewritten 

d 2 z d 

This corresponds to the motion of a particle in a potential V(,z) = e z . Using the initial condition 
([96]) which translates into (z — > — oo, dz/dt = 2) for t — > — oo, the first integral is 



This first order differential equation is readily integrated yielding 



(99) tanh" 1 W 1 - ^e z = t + C. 

Returning to original variables, we get 

where A is a constant of integration related to C. It is determined by ^(0) = yielding 8A 2 = 1 
so we finally obtain 

(101) e ~* = (TTW 

Using the Gauss theorem, the total mass is given by 

dc , , AM 

102 -it! = -— . 

dr 2ttR 

Let a = (XxPo/ Dy/ 2 R denote the value of £ at the edge of the box so that £ = ar/R. With 
these notations, the foregoing relation can be rewritten 

(103) anj/(a) = 

v ; ^ K J 2-kD 

Using Eq. (1 1001) and introducing the critical mass (|4T|) . we find that this relation is equivalent 
to 

a 2 M 



(104) 

v ' 8 M c -M 
Recalling the definition of a, this equation determines the relation between the mass and the 
central density according to 

( 105 ) A) = nR2 ,_M_ - 
Finally, noting that 

(106) p= $ , 

(i + f(£) 2 ) 2 

and using Eqs. (11041) and (11051) we finally obtain Eq. (jiljl . On the other hand, in an infinite 
domain, the Gauss theorem lim^ +00 £V>'(£) = and Eq. ffTOll imply that M = M c = Sgp. 
Then, using Eqs. (E3D and ffTUTT) we obtain Eq. (HBJ. 

The two-dimensional Boltzmann-Poisson equation (j57|) appeared in very different topics: 
self-gravitating isothermal gaseous cylinders in hydrostatic equilibrium [29], [30], statistical equi- 
librium states of two-dimensional stellar systems in the microcanonical ensemble [26j [31], sta- 
tistical equilibrium states of two-dimensional self-gravitating Brownian particles in the canon- 
ical ensemble [7], statistical mechanics of point vortices in two-dimensional hydrodynamics 
[4T1 [331 [341 153] . chemotaxis of bacterial populations [T8] . 



22 



D Explicit solution of the ID Boltzmann-Poisson equa- 
tion 



In this Appendix, we explicitly solve the Boltzmann-Poisson equation (I37p in d = 1. Writing 
p = po e ~^ with if) = /3m($ — $ ) where p is the central density and $ the central potential, 
and introducing the scaled distance £ = (2GPmpoY^ 2 x, the Boltzmann-Poisson equation in 
d — 1 can be written 

(-> 

with ^(O) = ^'(0) = 0. This is similar to the equation of motion of a particle in a potential 
V(i/j) = e~^. The first integral is E = ^(dip/dt;) 2 + e~^. Using the initial conditions we find 
that E — 1. Therefore, we get 

(108) ; ^ = dt 

>/2(l - e~+) 

This is integrated into 

(109) tanlT 1 y/l - = -^=£, 
and we finally obtain 

(110) 



cosh 2 ^/^)' 

In a infinite domain, the Gauss theorem (d<&/dx)(+oo) = GM in scaled variables becomes 

(111) lim ^'(0 = M(^) 1/2 . 

Using Eq. ffTTOl . this yields Eq. p7]) -b. Then, Eq. fTTTOl) yields Eq. (1391. In a finite domain, 
we call a = (2f3Gmpo) l l 2 R the value of the scaled distance £ at the box radius R. Then, we 
have £ = (a/R)x. The Gauss theorem (d$/dx)(R) = GM in scaled variables becomes 

(112) aip'(a) = pGMmR = rj. 
Using Eq. (IllOp . we find that 

(113) T) = V2at&nh(a/V2). 
We also have by definition 



;ii4) 



2Rp Q a 



2 



M 7] 

Therefore, Eqs. (11131) and (11 14ft determine the central density in terms of the temperature. 
Finally, the density and mass profiles can be written 

(115) p(x) = 5-^ — 7= — , M(x) = _ tanh(ax/V2R). 

V ' PV ' cosh 2 (ax/V2R) tanh(a/ v / 2) 

The central density monotonically increases as the temperature decreases (see Fig. H|). At 
T = 0, the density profile is a Dirac peak: p(x) = M5(x) (see Fig. [3]). Noting that a = 
PGmRa/ \/2, these results are equivalent to those of Sec. 15.21 The one- dimensional Boltzmann- 
Poisson equation (11071) appeared in very different topics: highly flattened galactic disks |48j, 
stellar systems stratified in plane parallel layers [4TJ], one dimensional self-gravitating Brownian 
particles [7], statistical mechanics of two-dimensional turbulence in a shear layer |50j . 
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